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ABSTRACT 



Using asteroseismic data and stellar evolution models we make the first detec- 
tion of a convective core in a Kepler field main-sequence star, putting a stringent 
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constraint on the total size of the mixed zone and showing that extra mixing 
beyond the formal convective boundary exists. In a slightly less massive target 
the presence of a convective core cannot be conclusively discarded, and thus its 
remaining main-sequence life time is uncertain. Our results reveal that best-fit 
models found solely by matching individual frequencies of oscillations corrected 
for surface effects do not always properly reproduce frequency combinations. 
Moreover, slightly different criteria to define what the best-fit model is can lead 
to solutions with similar global properties but very different interior structures. 
We argue that the use of frequency ratios is a more reliable way to obtain accu- 
rate stellar parameters, and show that our analysis in field main-sequence stars 
can yield an overall precision of 1.5%, 4%, and 10% in radius, mass and age, 
respectively. We compare our results with those obtained from global oscilla- 
tion properties, and discuss the possible sources of uncertainties in asteroseismic 
stellar modeling where further studies are still needed. 



Subject headings: asteroseismology 
oscillations — stars: interiors 



stars: fundamental parameters — stars: 



1. Introduction 



A detailed comprehension of the physical processes taking place in deep stellar inte- 
riors is of paramo unt importance for an accurate description of stellar populations (e.g., 
Chiosi et al.l Il992l ). Assumptions and parameterizations used to represent processes like 
convection, rotation, overshooting and microscopic diffusion directly impact the quantities 
normally used to characterize different types of stars, such as effective temperatures, col- 



ors, surface gravities and composition (see e.g.. lRoxburghlll978l IVandenberg fc Smithlll988 



Michaud &: Proffittlll993l . lMaeder fc Meynetll200Cll . and references therein). As a consequence 
it is extremely difficult to determine reliabl e stellar ages of fi eld main-sequence stars, where 
uncertainties are well above the 20% level ( ISoderblomI 120101 ). 



One aspect of critical relevance for dating main seq uence stars is the existence and 
size of a convective core. It has been long known (e.g., iMaeden Il974j ) that reproducing 
the observed hook-like feature in the Color- Magnitude Diagram (CMD) of intermediate- 
age clusters requires that stars at the turn-off phase possessed a convective core during 
their main-sequence lifetime. Since the metallicity of cluster stars is assumed to be fixed 
and homogeneous, the age of the clust er can be estimated by determining the bes t fit of 



isochrones to the observed CMD (e.g.. IVandenberg &: Belli Il985l . ISalaris et al.l 119971 ). The 



great caveat of this technique is that the critical mass at which a convective core appears 
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in evoluti onary models, normally around 1.1 Mp ^ , strongly depend s on the input physics 



[see e.g., iChristensen-Dalsgaard &: Houdekl l201Cll . iMagic et al.l I2OIOI ). Moreover, the total 
size of the mixed central region in stellar models can be controlled by a parameter of mixing 
efficiency beyond the formal convective boundaries of the core; this is calibrated to match 
observations of clusters a nd is again dependent on the rest of the input physics used to 
construct isochrones (e.g.. iMaeder fc Meynetlll99ll . iPietrinferni et al.ll2004j ). Therefore, we 
currently lack a reliable estimation of the size of convective cores in main-sequence stars. 

An exciting approach to overcome these impediments comes from the possibility of 
piercing the surface of stars by studying their pulsations, and the exquisite data currently 
being obtained by the Kepler mission provide an excellent opportunity to achieve this 
goal (IChristensen-Dalsgaard et al.l 120071 iBorucki et al.l l2009l ) . Several theoretical studies 
of how asteroseismology can yield information about stellar cores have been carried out (e.g.. 



2007 



Audard &: Provost 1994. Mazumdar fc Antia 



Cunha fc Brandao 



2011 



2001 



Silva Aguirre et al. 



Mazumdar et al.ll2006, Cunha fc Metcalfe 



2011aJ), but the paucity of data for main- 



sequence targets prior to the launch of Kepler allowed very few investigations of this type in 
stars other than the Sun. In fact, all previous da t a-based asteroseismic s t udies of stellar cores 



prop erties (such as iMiglio fc MontalbanI |2005| . iDeheuvels et al.1 l2010l . Ide Meulenaer et al. 



2OIOI ) relied on ground-based observations of much shorter time coverage and lower qual- 
ity than the multi-month Kepler observations. 

From the hun dreds of main-sequen ce targets for which the Kepler mission has provided 
asteroseismic data (jChaplin et al.ll201ll ). individual frequencies ha ve now been determined fo r 
dozens of them from time-series spanning several months (e.g., lAppourchaux et al.ll2012l ). 
We carry out here the first investigation using Kepler data on two of these targets with 
the aim of detecting convective cores, constraining their size, determining accurate stellar 
parameters, and narrowing down the age uncertainties compared to those obtained by fitting 
stellar tracks. 

To accomplish this, we proceed as follows: in Section |2] we describe the available seismic 
and spectroscopic observations, which are used in Section|3]to obtain initial stellar parameters 
from global asteroseismic fitting techniques. The large uncertainties obtained from this 
analysis, particularly in age, lead us to consider diagnostic tools sensitive to the deep interior 
of stars, presented in Section HI Several teams performed detailed modeling of the targets 
using different evolution and pulsation codes, described in Section |5l The results of this 
detailed analysis are given in Section |6l together with the derived stellar parameters. We 
discuss our findings, as well as other possible sources of uncertainties, and give closing 
remarks in Section |71 
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2. Target selection and observations 



Table 1: Parameters determined from the observations of the Kepler mission, ground based 
photometry, and spectroscopy for both stars. Luminosity of Perky determined from par- 
allax measurement s and bolometric corrections. Kepler magnitude Kp as determined by 
Brown et all J201l[ ). 



Star 



KIC 



Kt, 



I'max (/iHz) (Au) (flUz) T^ff (K) 



log 9 



[Fe/H] 



log(L/L0) 



Perky 
Dushera 



6106415 
12009504 



7.18 
9.32 



2210 ± 50 
1833 ± 40 



104 ±0.5 
88 ±0.6 



6000 ± 200 
6200 ± 200 



4.27 ±0.1 
4.30 ±0.2 



-0.09 ±0.1 
0.0 ±0.15 



0.26 ±0.04 



Based on the full ensemble of stars with 9 months of Kepler observations available, we 
initially selected targets with seism ic characteristi cs belonging to the main-sequence phase 
(no signature of mixed- modes, see iBeddingI l201ll . for a thorough explanation). Since we 
are interested in detecting convective cores, we further discarded the targets with global 
seismic properties suggesting a mass value lower than solar, and finally chose two stars 
from the r emaining targets with a large number of reliable individual frequencies deter- 
mined (see lAppourchaux et al.ll2012l ). These are dubbed Perky (KIC 6106415) and Dushera 
(KIC 12009504). 

The power spectrum of solar-like oscillators is modulated b y a Gaussian-like env elope, 
and presents a near-regular pattern of overtone frequencies (e.g., IChaplin et al.ll201ll ). Two 
parameters can be readily extracted from it, namely the frequency of maximum power of 
oscillation z/max and (Au) , the average separation between the consecutive p-mode overtones: 

Aue{n) = Un^f, - Un-i^e , (1) 

where Un/ is the mode frequency of angular degree i and radial order n. These are usually 
referred to as the global oscillation parameters. 



Values of (Au) and i^max were obtained for both tar gets as described by 



f|2011l) using time series prepared from short-cadence data (JGilliland et al. 



2010 



Verner et al. 



Garcia et al. 



20 111) . Individual frequencies were extracted in the manner explained by lAppourchaux et al. 
(120121 ). resulting in 33 oscill ation modes for Perky and 34 for Dushera. These are listed in 
Table [7] and in Table A. 53 of lAppourchaux et al.l (120121 ) for Perky and Dushera, respectively. 



Effective temperature T^s, surface gravity log g and iron abundance [Fe/H] were avail- 
able for b oth targets from hig h-resolution spectroscopic observations analyzed with the VWA 
software (JBruntt et al.ll2012l . and references therein). The case of Dushera presented a chal- 
lenge: a set of low-resolution spectra was analyzed with a standa rd method based on the 
equivalent widths of iron lines using ARES and MOOG (see e.g., ISousa et al.ll2007l ). This 
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method, described in detail in ISousa et al.l (|2008[ ). yielded atmospheric parameter values 
that were not consistent with those obtained from VWA. Although one might be tempted to 
disregard the results from the l ower quality data, estirnations of Tpfp for Dushera using the 
InfraRed Flux Method (IRFM. ICasagrande et al.ll2010l . ISilva Aguirre et al.ll2012l ) suggest a 
value that is consistent with that obtained by ARES+MOOG. Since a thorough comparison 
of these spectroscopic methods is beyond the scope of this paper, we adopted average values 
and large uncertainties in the observed parameters that encompassed both determinations. 



For the case of Perky only the results from iBruntt et al.l (120121 ) were available and 
the IRFM Teff determination is consistent with this value. However, based on our experi- 
en ce with Dushera, we decided to use more conservative uncertainties than those quoted 
by iBruntt et al.l (120121 ). Since parallax measurem ents are also availa ble for this star, we 
estimate d its luminos i ty usiiig the Hipparcos data (jVan LeeuwenI l2007l ) and bolometric cor- 
rections (lFlowerlll996l . lTorresll2010l ). The complete set of observational constraints considered 
in this study is given in Table [H 

An initial determination of ages can be ma de using these spectr oscopic parameters and 
stellar isochrones following a bayesian analysis (ISerenelli et al.ll2013l ). For the values quoted 
in Table [H this technique yields age uncertainties of 50% and 60% for Perky and Dushera, 
respectively. Even when reducing the errors in Teg to 100 K and considering and error in 
log g of 0.1 dex for both targets, the obtained uncertainties only decrease to 36% and 50%. 
Ages determined with a precision of ~20% can only be achieved with this technique for 
targets with accurately determined observational constraints, w 'hich are not availab le for 
these relatively faint field main-sequence stars (see Section 5.2 in lSerenelli et al.l ( l2013l )). In 
fact, age results for main-sequence intermediate- mass field stars have often uncertainties of 
20-50% when determined from a single ste llar evolution data base, that is, before systematics 
are taken into account (see Section 4.1 in ISoderblomI (120101 )). 



3. Global asteroseismic analysis 

Using the parameters given in Table [H it is possible to apply global asteroseismic tech- 
niques and determine the stellar prope rties of the ta rgets. (Au) approximately scales as the 
square root of the mean density (e.g., IUlrichlll986l ). while z/max seems to be related to the 
acoustic cutoff frequency of the atmosphere and thus to the surface grav i ty and effective 



temperature of the star (e.g., iBrown et al.lll99ll . iKjeldsen fc Bedding! Il995l . iBelkacem et al. 



201ll ). Based on these dependencies, two scaling relations anchored at the solar parameters 
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Table 2: Stellar parameters as determined from each pipeline for the two targets. Also shown 
are the values obtained with the direct method. See text for details. 



Sta 



Perky 



Parameter 



Marigo 



YREC 



RADIUS 



GARSTEC 



RadExlO 



M(Mq) 

R(Rq) 

log S 

Age (Gyr) 



1.06 
1.21 
4.30 
5.73' 



+0.09 

0.07 
+0.03 

0.03 
+ 0.01 

0.01 
+ 2.38 

2.44 



1.06 
1.21 
4.30 
6.00 



+ 0.0S 

0.07 
+ 0.03 

0.03 
+ 0.01 

0.01 
+ 2.44 

2.18 



1.12 
1.23 

4.31 
3.71 



+0.05 
-0.06 
-0.02 
-0.02 
+ 0.01 
0.02 
+ 1.80 
1.22 



1.09 
1.22 
4.31 
6.21 



+0.07 
-0.06 
-0.03 
-0.02 
-0.01 
-0.01 
+ 2.24 
-2.15 



1.10 ± 0.08 
1.22 ± 0.03 
4.30 ± 0.01 
6.77 ± 2.20 



1.24 
4.30 



,+0.10 
^-0.10 
-0.01 
-0.04 
+ 0.01 
0.01 



1.03' 
1.20 
4.29 



+0.10 
-0.10 
-0.04 
-0.03 

+ 0.02 
0.02 



1.08 ± 0.09 
1.22 ± 0.03 
4.30 ± 0.01 



4 21+3. 3S 
-2.21 



5.52+^;^* 6.37 + 3.18 



1.07 
1.22 
4.30 
5.80 



+ 0.06 

0.10 
+ 0.02 

0.04 
+ 0.01 

0.01 
+ 3.70 

2.20 



Direct 



1.10 ± 0.07 
1.23 ± 0.02 
4.30 ± 0.01 



1.25 
1.43 



+ 0.08 

-0.10 

+ 0.04 

0.04 



1.21 



+ 0.09 

0.08 

, + 0.03 



Dushera M(Mq) 

R{Rq) 
log 9 
Age (Gyr) 3.2lt\-ll '■i.82t\ll 



1.23 



+ 0.07 
0.05 



1.27 
1.44' 



+0.10 

0.10 

+ 0.04 

-0.04 



4 22+0-01 4 22+001 
^■^^-0.01 ^■''^-o.oi 



4 22+0-01 
^■^^-0.01 

2 94+1-06 
^■^^-1.05 



4 22+0-01 
*-^^-o.oi 

3.78tlil 



1.21 ± 0.10 
1.41 ± 0.05 

4.22 ± 0.01 
6.43 ± 2.06 



,.+0.09 
"-0.09 



4 22+0-01 
^■^^-0.01 

2.89i;;0i 



1.16 
1.41 



+ 0.14 

-0.13 

+ 0.05 

0.05 



1.26 + 0.09 1.19Io.o8 1.22 + 0.08 
1.44 + 0.04 1.41+0-0^ 1.42 + 0.03 



4.21+0-0^ 4.22 + 0.01 
^■''^-I'll 3.74+1.71 



4.22^0.0^ 4.22 + 0.01 

4 05+l-**l 
*-O'-'-1.05 



can be written 



M 
R 

Rq 



iy„ 



'max,0 



'^max,© 






-4 



-2 



T, 



cflf 



^eff,0 
^cff 



z 



cff,0 



3/2 



1/2 



(2) 



(3) 



Here, T^s^q = 5777 K, At/^^ = 135.1 ± 0.1 /xHz and z/, 



inax,0 



3090 ± 30 yuHz are the observed 
values in the Sun (JHuber et al.ll201ll ). Provided a measurement of Tgff is available, these 
relations give a determination of mass and radius independent of stellar models (the direct 
method, see e.g.. lMiglio et al.ll2009l . IChaplin et al.ll201lUSilva Aguirre et al.ll2011bf ). It is also 
possible to search for a best fit to these parameters within pre-calculated sets of evolutionary 
tracks, where chemical composition can be taken into account (the gr id-based method, e.g., 
Stello et allbood lOuirion et allboiol . ICai et al.ll201l[ ISasu et al.ll20li ) . The latter approach 
also provides a determination of the stellar age. 

We performed a grid-based search using several available sets of evolutionary models 
to make an initial estimation of the stellar parameters of both targets. Details on how the 
sets of evolutionary tracks w ere construc t ed an d the different s e arch methods can be found 
in the following references: ISasu et al.l (l2010h and iGai et al.l (1201 ih for the YY, YR EC, 
Dotter and Marigo results: Istello et al.N2009[ ) for R ADIUS: Isilva Aguirre et al.l J2012I ) for 
GARSTEC; ICreevev et all (]2012|) for RadExlO; and lOuirion et al.l (I2OI0I) for SEEK. One 



grid-based pipeline built using the BaSTl|^ isochrones ( Pietrinferni et al.ll2004J ) is presented 
here for the first time. 



Briefly, the BaSTI pipeline uses a denser set of isochrones in metallicity space than the 
publicly available ones, computed specially for the latest revision of the Geneva-Copenhagen 



^http://albione.oa-teranio. inaf.it/ 



Survey flCasagrande et al.M2011l) and n ow ex tended for asteroseismic analysis. They are 
built assuming the iGrevesse fc Noelsl (Il993l ) solar composition and adopting a value of 
Ay/AZ = 1.45, in agreement with Big Bang nucleosynthesi s (SBBN) values of p rimordial 
hehum abundance at low metalhcities (Yp = 0.2482 ± 0.0007. IsteigmanlboOTl . I2OI0I ). The 22 
available chemical mixtures cover a metallicity range between —3.27 <[Fe/H]< +0.5. Models 
used here do not include ov ershootirig and gravitational settling, and are calculated assuming 
mass loss according to the iReimersi (119751 ) formulation with a fixed efficiency parameter of 
1] = 0.4. The algorithm to find the best match to the data is that of the GARSTEC grid, 
based on Monte Carlo realizations to form the probability distribution of each parameter. 

Table [2] shows the results obtained with each grid pipeline for the two targets, including 
their 1-a uncertainties. As mentioned above, some methods return a probability distribution 
for the stellar parameters characterized by asymmetric error bars. The last column shows, 
for comparison, the values obtained with the direct method (Eqs. [2] and [3]). There is an 
overall good level of agreement in the results for log g and R/Rq, explained by the fact that 
the global asteroseismic parameters are mostly sensitive to ratios between stellar mass and 
radius (mean density and surface gravity). Due to these st rong dependen c ies, stellar radii 
and l og g as determined from seismology are very precise (JBeddingj I2OIII iMorel fc Miglio 
2012l . lSilva Aguirre et alll2012[ Huber et al. 2012). 



On the other hand, asteroseismic mass determinations still require verification from 
independent methods. The mean values from the grids are 1.08 Mq for Perky and 1.23 Mq 
for Dushera. All masses obtained by the 9 pipelines are contained within a 10% range from 
these averages. In spite of that, when taking into account the 1-a uncertainties given by 
each grid, mass values between ~0. 95-1. 20 Mq and ~1. 05-1. 35 Mq are possible for Perky 
and Dushera, respectively. The large age differences in the results (~50% or more) are a 
direct consequence of this mass difference (~0.3 Mq). Moreover, the fact that the results 
are in the vicinity of the ~1.1 Mq limit where the onset of core convection is thought to 
occur suggests that the existence and extent of a convective core will play a significant role 
in the internal structure of the star. 

These results are to be taken with some caution since the scaling relations given in 
Eqs. [2] and [3] are not exact. Moreover, they are of course dependent on the uncertainties in 
asteroseismic ((Az/) and t'max) as well as atmospheric (Tes and metallicity) inputs fed into 
the grid searches, which we decided to keep at a conservative level as explained in Section [2l 
In fact, each grid of models has been constructed using different sets of input physics, solar 
mixture, opacities, equation of state, etc., which accounts for some of the differences seen in 
the results (see Chaplin et al. (2013), in preparation, for a throughout discussion of the grid- 
to-grid scatter). When bearing this in mind it is clear that a detailed asteroseismic study 
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is mandatory to determine masses with higher accuracy, to explore the energy transport 
mechanisms in the inner structure, and to better constrain the age of the targets. The lack 
of sensitivity of the global seismic analysis to the finer details of the stellar structure needs 
to be compensated with tools capable of probing the deep interior of stars. 



4. Asteroseismic diagnostics for stellar interiors 

Different combinations of oscillation frequencies have been proposed in the literature to 
study the internal structure of stars. The two most commonly used are the large frequency 
separation (Eq. [T]) and the small frequency separation defined as 



di^i+2{n) = Vn,(. - V, 



n-l/+2 



(4) 



The latter quantity is sensitive to the conditions in the stellar core and thus to the cen- 
tral hydrogen content (age) for main -sequence stars (e.g.. IUlrichlll986l . lGoughlll987l . iGabriel 
19891 . [Roxburgh fc Vorontsovlll994bl ). In contrast, the l arge frequency separation is known 



to be affected by the outer layers of stars (see e.g.. IChristensen-Da 



Christensen-Dalsgaard &: Perez Hernandezlll992l . [Roxburgh &: Vorontsov 



sgaard et al 



2003 



1988 



Roxburghll2010f ). 



where stellar models use simple versions of convectio n theory instead o f properly dealing with 



turbulent pressure and non-adia batic effects (e.g., iBalmforthl Il992l . [Rosenthal et al.l Il999 
Ballot et al.ll2004j . lHoudekll2010[ ). For the solar case, differences in Au around the y ma;x.R, 
region between Bi SON (Birmingham Solar Osc i llatio ns Network, IChaplin et al.l 119991 ) data 
and Model S from lChristensen-Dalsgaard et al.l (119961 ) are of the order of 1/iHz, much larger 
than the typical solar observational uncertainties of a^^, ~0.06/iHz. 

In order to study the deep interior of the star we need combinations of frequencies that 
isolate the signal arising from the stel lar center as much as possible from surface contami- 
nation. iRoxburgh fc Vorontsovl ( 120031 ) proposed to use ratios of small to large separations . 



which are mainly deterna i ned by the inner structure of the star ( IRoxburgh fc Vorontsovll2003 



Oti Floranes et al.ll2005l . lRoxburghll2009l ). They are constructed as 

. . do,2{n) 



roi(n) 



ro2[n) 



doi{n) 



Az/i(n) 



'^lO 



dio{n) 



Ai/i(n)' ' "" Az/o(n + l)' 
where doi{n) and dio{n) are the smooth 5-points small frequency separations: 



(5) 



(6) 



doiin) = -(z/„-i,o - 4i/„_i,i + 6i/„,o - 4z/„,i + z/„+i,o) 



(7) 
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dlo{n) = --{jyn-1,1 - 4:iyn,0 + ^Vn,l - 4z/n+l,0 + T^n^l,l) 



(8) 



The ratios of degree and 1 (Eq. pi) have already been used to identify the location 
of the conve ctive envelope in the Sun ( Roxburgh! 200S), an d also in other main- sequence 
targets (e.g-. lLebreton fc Goupilll2012l . iMazumdar et al.ll2012h . Inspection of Eqs. [THH] shows 
that roi(n) and rio(n) are alternatively centered in Vn^ and z/„,i, respectively. Therefore we 
can consider the ratios roi(ra) and rio(n) as one unique set of observables as a function of 
frequency which we call roio: 



'"oio = {?^oi(^),^io('^),'"oi('^ + l),no(^ + l),roi(n + 2),rio(ra + 2),...} 



(9) 



Existence of a growing convective core in stellar models produces an abrupt discontinuity 
in the chemical composition at the boundary between the core and the radiative envelope. 
The formal location of this boundary is determined by a stability criterion, most commonly 
the Schwarzschild criterion. When additional mixing such as overshooting is included in 
stellar models, the transition between the homogeneously mixed center and the radiative 
zone does not occur at the formal edge of the convective core but at a certain distance from 
it determined by the adopted mixing efficiency. In either case, the discontinuity in chemical 
composition (and thus adiabatic sound speed) produces a shift in the frequencies relative to 
those in a corresponding star with a smooth sound-speed profile that is an oscillatory function 
of the frequency itself, and whose period is related to the loca tion of the discontinuity (e.g., 
Roxburgh fc Vorontsovl[i994al . l2nn7l . Icunha fc Metcalfelbon?! ). 



It is customary to represent the radial coordinate in stellar interiors by the wave travel 
time from the surface to the center, called the acoustic depth r, or its alternative represen- 
tation named the acoustic radius: 



t 



'^dr 

Cs 



(10) 



where Cg is the adiabatic sound speed. This quantity measures the travel time from the 
center towards the surface. Considering that the total acoustic radius is given by ttot = 
^(-Rtot), then r = ttot — i- If the sound-speed discontinuity is located at, say, ri in radial 
coordinates, and that same position is represented in acoustic radius and acoustic depth 
by ti and ti respectively, the periods of the oscillation in frequency space induced by this 
sharp variation are l/(2t i) and l/(2ri), according to the considered seismic coordinate (e.g., 
Montgomery et al.l 120031 ). 



We stress the fact that the frequency combinations will be sensitive to the position of the 
transition between the homogeneously mixed and radiative region, regardless if this coincides 
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or not with the formal boundary of the convective core. If the ratios were only affected by 
the location of the outer edge of the mixed region, information about its location could 
be directly extracted from the asteroseismic data. Nonetheless, the small acoustic radius 
characterizing the size of the convective core introduces an oscillatory signal with a period 
much larger than the observed frequency range. Moreover it has been shown that, although 
the ratios are mostly sensitive to the pres ence and size of a convective core, they are also 



affected by the central hydrogen content (JBrandao et al.l l2010l . ISilva Aguirre et al. 



2011a 



and still contain the signal of the base of the convective envelope (JRoxburgh fc Vorontsov 
20031 . lRoxbu7ghll2010h . 



There are several other ingredients of stellar evolution that play more or less important 
roles in the determination of the sound speed. Use of different opacities, equation of state, 
and convective theory applied affects the internal temperature stratification in stellar models. 
When using frequency ratios one must bear in mind that their behavior reflects the interplay 
between evolutionary stage (through the central hydrogen content) and a central convective 
region (when it is present), while the rest of the microphysics affects them to a much lesser 
(but not negligible and still not fully quantified) extent. A detailed study of these issues goes 
beyond the scope of this paper and will be presented elsewhere (Silva Aguirre et al. 2013, 
in preparation). 

Figure [1] shows the observed large frequency separations and ratios roio for the two tar- 
gets considered in this study. Inspection of the lower panel shows that the stars have different 
mean values and slopes in their ratios. Moreover, the oscillatory component introduced by 
the position of the convective envelope is clearly visible in the ratios of Perky. 

The behavior of the frequency ratios in both targets shows a sudden increase at high 
frequencies, that occurs regardless of the length of the data used and is not predicted by 
models. We have compared individual frequencies extracted from datasets of different lengths 
and confirmed that, due to the low signal-to-noise ratio and the larger linewidth at the high- 
frequency end, this behavior is the result of realization noise. For this reason, we will not 
consider the frequency ratios on this regime when comparing to stellar models. In the lower 
panel of Fig. [T] we m ark the frequency range th at will be taken into account to study tqio 
(the linear range, see ISilva Aguirre et al.l 1201 lal ). The range considered to fit the ratios ro2 
has also been restricted following the same criterion. 

Together with (Az/) the average small frequency separation is used to construct a n 
asteroseismic HR diagram known as the C-D diagram ( IChristensen-Dalsgaardl Il984l Il988l ). 
A slightly modified version of this diagram is constructed by using the ratios ro2 instead 
of the small frequ ency separations, because it is a more effective indicator of stellar age 
( jWhite et al.ll201l[ ). In Fig. [2] we compare the position of Perky and Dushera in this diagram 
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with evolutionary tracks at solar metalliticy, showing that for a given composition and input 
physics the age of the targets is well constrained by ro2. Also plotted in the figure are tracks 
with a metallicity ~0.2 dex lower, showing that the offset due to composition has a strong 
impact in the derived ages from this diagram. This is particularly important for Dushera, as 
its position is closer to where the tracks collapse in this diagram when entering the subgiant 
phase of evolution. 



110 




1400 1600 1800 2000 2200 2400 2600 2800 

u (/J.Hz) 



Fig. 1. — Upper panel: large frequency separations Ai/q, Ai/i, Az/2 of Perky (filled red di- 
amonds) and Dushera (open blue circles). Lower panel: frequency ratios roio('^)- The 
frequency range selected to fit the ratios is depicted by horizontal dotted and dashed lines 
for Perky and Dushera, respectively. See text for details. 
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Fig. 2. — Modified C-D diagram showing the location of Perky (filled red diamond) and 
Dushera (open blue circle). Evolutionary tracks are shown (from right to left) for masses 1.0- 
1.3 Mq at solar [Z = 0.017, solid lines) and sub-solar {Z = 0.011, dashed lines) metallicities. 
Dotted line s connect pos i tions at the same age for the solar metallicity tracks. Evolutionary 
tracks from Iwhite et al.l (1201 ll ) . 



5. Detailed modeling 

Asteroseismic modeling of individual targets has experienced a lea p forward thanks to 



the accurate frequencies extracted from CoRoT and Kepler daia. (e.g,., iBarban et al.l 12009 



Ballot et al.l I20 111 ICampante et al.l 1201 ll . lAppourchaux et al.ll2012l ). The usual strategy ap- 
plied in these type of studies consist on calculating different sets of one-dimensional hydro- 
static stellar models. Using these as the equilibrium structures for a code capable of comput- 
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ing frequencies of oscillations, it is possible to produce a theoretical set of the spectroscopic 
and asteroseismic observables. Finally a criterion mus t be applied that defines which set 
of ca l culations best reproduces the available data (e.g. 
201ll . iMathur et aPboid IPoean et allboiSJ ). 



Metcalfe et al.ll20ld . iBrandao et al. 



There are many ingredients in the computations of stellar models that are still not fully 
constrained, both in the microphysics (e.g., equation of state, opacities, nuclear reactions) 
and in the macrophysics (e.g., conve ctive energy transport) of stellar evolution (see e.g., 
Cassisil l2009l iKippenhahn et al.l |2013| ) . Even when the same assumptions are used in dif- 
ferent evolutionary codes, the resulting niodels can differ in the ir global properties, internal 
structure (see e.g 



Le 



jreton et a 



frequencies obtained ( iMoya et al. 



20081). 



20081 . iMontalban et al.l |2008[ ) , and also in the theoretical 



For these reasons, several teams performed detailed modeling to match the stellar pa- 
rameters presented in Table [1] and individual frequencies, using the stellar evolution and 
pulsation code of their choice (see Section [STTj below) . This allows us to obtain an at least 
qualitative measure of the systematic uncertainties inherent in the full modeling process, 
arising from the different input physics and implementations of search methodologies. A 
quantitative investigation of the individual influence of the various differences between the 
models presented here is beyond the scope of this paper. 

Techniques to find the b est-fit model to a given set of observables usually include mini- 



miza tion of a reduced x^ (e.g.. lDeheuvels &: Michelll201ll . lBrandao et al.ll201lllGruberbauer et al. 



20121) defined as 



X 



N ^ 






(11) 



where k runs over the observable parameters, Xobs(^) represent each observable with its 
uncertainty cr(a;obs(^)), a;niodei(^) are the theoretically produced values, and A^ the total 
number of observations. In the analysis of the models that best fit the data, we will consider 
X^ values for individual frequencies and ratios. It is clear from Eqs. [7] and [8] that the 
differences entering the ratios roio are strongly correlated. We take this into account by 
including their covariance term in the x^ 



X 



^ l^^obs 



•^modclj 



\T 



l^'^'obs •''model y 



(12) 



where C is the covariance matrix calculated after perturbing the observed frequencies using 
Monte Carlo simulations. To calculate the goodness of fit for the ratios tqio we use Eq. [TH 
which reduces to Eq. [TT]for the individual frequencies and ratios ro2 as they are independent. 

As mentioned in Section |U theoretical frequencies of oscillations are not able properly to 
reproduce the observed ones due to our poor understanding of the physical processes taking 
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place in the outer layers of stars (see e.g., IChristensen-Dalsgaard et al.lll996l . for the solar 
case). There are different formu lations to correct for this offset (see e.g., the supplementary 
materials forlCarter et al.ll2012l . and references therein), the most commonly used is that of 
Kjeldsen et al.l ( 120081 ). It consists of an empirical power law derived using helioseismic data 
and standard solar models, which can be written as 



z^obs(n) -rure{{n) = a 



^0 



(13) 



Here z/obs('^) correspond to the observed frequencies and z/ref(Ti) are the theoretical ones. To 
apply this surface correction one has initially to specify the values of the exponent b and the 
reference frequency uq. Then r is determined after calculating an observed and theoretical 
average of the large frequency separation, obtained using a linear fit to the individual i = 
modes as a function of radial order n. The value of a is determi ned from the averag e 
theoretical and observed i = frequencies (see Eqs. 6, 8, 9 and 10 in lKjeldsen et al.ll2008l ). 
Although most modeling teams use this formulation in their analysis, they determine the 
parameters involved in the correction in slightly different manners (see Section \572\ below) . 

Another point where discrepancies can arise from is the estimation of the surface iron 
abundance in stellar models, to be compared with the [Fe/H] value determined from spec- 
troscopic observatio ns. Although the issue of the solar surfa ce abundances is still a matter 



of debate (see e.g., iGrevesse et al.l l2010l . ICaffau et al.ll2010f ). the observed [Fe/H] in stars 



barely depends on the considered solar reference because the largest discrep ancies between 



the d ifferent solar mixtures come from the C, N, O and Ne elements (e.g., ISerenelli et al. 



20091 ). On the other hand in stellar models it is assumed that 



\og{Z/X) - logiZ/X)^ ^ [Fe/H] , 



(14) 



where X and Z are respectively the mass fraction of hydrogen and all elements heavier 
than helium, and log{Z/X)Q is the surface solar value. Equation [H] clearly depends on 
the chosen solar mixture, and assumes that the fraction of each element comprising Z is 
distributed according to the selected set of solar abundance ra tios. The commonly ado pted 
values of {Z/X)p-, vary for different mixtures between 0.0245 (IGrevesse fc Noelslll993l ) and 
0.0181 (jAsplund et al.ll2009l ). This means that for models with the same (Z/X) at the surface 
but computed with a different set of solar abundances, the theoretically determined [Fe/H] 
can vary by ~0.13 dex, comparable to the uncertainties quoted in Table [T]from spectroscopic 
determinations. 

Due to the issues aforementioned, and the large error bars in the spectroscopically deter- 
mined [Fe/H] and T^s (see Tabled]), we will not include these constraints in the calculations 
of goodness of fit via x^- These observables were used by all teams as guidelines when 
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exploring the parameter space (except for models G, see Section [5.1.71 below), and will be 
considered as reference for comparison that only in extreme cases of discrepancy could be 
used to discard models. 



5.1. Evolutionary calculations and search algorithms 

In the following sections we describe the evolutionary code, input physics, and fitting 
techniques adopted by each modeling team to find the best match for both targets. 



5.1.1. Models A 
Models A were c omputed using the Garching Stellar Evolution Code (GARSTEC, 



Weiss fc Sch 



Rogers et al 



1996 



attlll2008[). We adopted the 200 5 version of the OPAL equation of state (EOS, 



Rogers &: Navfonovl 120021) c omplemented by the MHD equa tion of state 



for low temperatures ( Hummer fc Mihalaslll988l ). l ow-temperature opaciti es fro mlFerguson et al 



(2005) and OPAL opacities for high temp eratures (jlglesias fc RogersI 



19961 ). the lGrevesse fc Sauval 



( Il998l ) solar mixture, and the NACRE (jAngulo et al.l 119991 ) compil ation for thermo r iuclea r 
reaction rates with the updated cross section for ^^N(p, 7)^^0 from iFormicola et al.l ( l2004l ). 
Due to the large uncertainty in the surface composition of Dushera (see Table [1]), diffu- 
sion of heliu r n and heavy elements was included only for the models of Perky using the 
Thoul et al.l (119941 ) prescription. The atmospheric stratification is given by Eddington's 
gray plane-parallel formalism. 

Convective zo r ies ar e treated with the mixing-length theory (MLT) as described by 
Kippenhahn et al.l (|2013[ ). using the solar calibrated value of the convective efficiency Omit = 
1.791. Overshooting is implemented in GARSTEC as a diffusive proces s consisting of an ex - 
ponential decline of the convective velocities within the radiative zone (JFreytag et al.lll996l ). 
The diffusion constant is given by 



Dov [z) = Do exp 



-2z 



(15) 



where ^ corresponds to an efficiency parameter, Hp is the pressure scale height, z the distance 
from the convective border, and the diffusion constant Dq is derived from MLT-convective 
velocities. Using open clusters, the calibrated value for the overshooting efficiency is defined 
as .^ = 0.016. To prevent the survival of the pre-main sequence convective core due to 
overshooting in the solar model, a geometric factor for efficiency restriction is used when 
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convective cores are small (jSchlattl fc WeissI Il999l . iMagic et al.ll2010l ). When a convective 
core was present, we explored overshooting efficiencies between ^ = 0.004 — 0.035. 

Oscillation frequencies of the specific models analyzed were compute d using the Aarhus 
Adiabatic Oscillation Package (ADIPLS, IChristensen-DalsgaardI l2008al ). Models were se- 
lected by matching the observed large frequency separation and minimizing the x^ fit to 
both frequency ratios roio and ro2. No surface correction was applied to define the best fit 
to the data. 



5.1.2. Models B 



Models B were computed using the Aarhus Stellar Evolution Code (ASTEC. IChristensen-Dalsgaard 



2008bl ) . The models are calculat ed using the NACRE compilation of nuclear reaction rates 
and the iGrevesse fc Noeld (119931 ) solar mixture. The EOS is the OPAL 2005 versio n, while 
the opacities are those of OPAL for high temperatures and iFerguson et al. ( 20051) for low 
temperatures. Convection is treated using the MLT as described by iBohm-Vitensd (Il958l ). 
and no microscopic diffusion is considered. 

Core overshooting is implemented as an extension of the central mixed region using the 
radiative temperature gradient, and controlled by an efficiency parameter Oov 



cin 



aovmin(rcc, Hp) 



(16) 



where Tr. 



is the radial size of the convective core. This prescription restricts overshooting 
if convective cores are smaller than Hp. Values of ctov = 0.0,0.1,0.2 were explored in the 
calculations of models B, similar to the calibrated efficiencies of other codes with the same 
prescription (see models F below). 



The atmospheric stratification is given by the T— Top-relation of model C in lVernazza et al 



(Il98ll ). Oscillation frequencies are computed using ADIPLS to find the best ffi to the data 
by minimizing the y^ va lues of individual frequencies and frequency ratios after applying 
the iKjeldsen et al.l (120081 ) surface correction. There is one subtle difference in the way the 
surface correction is applied in models B, that will turn out to be of particular importance 
in the case of Perky (cf. Eq. [T^ . To obtain the average large frequency separation needed 
to compute r, a ffi is made to all available frequencies instead of only the radial modes. 
Similarly, radial and non-radial modes are used in the average frequencies used to obtain the 
a parameter. 
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5.1.3. Models C 



Models C were also computed using ASTEC with the same input physics as described 
for models B, but no overshooting was included. A grid of models was calculated covering the 
following range in the mass-initial composition space: M= 1.0 — 1.6 Mq, Y = 0.24 — 0.32, and 
Z/X = 0.01 — 0.07, while the mixing length parameter omlt "was fixed to 1.8. Frequencies 
were computed with ADIPLS for the models that have [Fe/H] and log g within 3-0", and T^s 
within 1-0" from the values given in Table [H since the uncertainties on Tes were relatively 
conservative. The best-fitting models were selected based on minimizatio n of the x^ value 
for only the individual frequencies after applying the iKjeldsen et al.l (120081 ) empirical surface 
correction. 



5.1.4. Models D 



Models D were calculated using the CESAM2K evolutionary code (Morel 19971 . iMorel fc Lebreton 



20081), with the OPAL 2005 EOS a nd opacities, the 



volutionary code (Morel 19971 . iMorel fc 
Grevesse fc Noeld (119931 ) solar mixture 



and mi croscopic diffu s ion fo llowing iMichaud fc Proffittl (119931 ). Convection is modeled us- 
ing the ICanuto et al.l (119961 ) formulation, which for this case gives a solar calibrated value 
for the convective efficiency of aconv = 0.64. Nuclear reactions are those from the NACRE 
compilation, while the atmosphere is described by Eddington's gray law and connected to 
the stellar interior at an optical depth of Top = 10. Core overshooting is implemented using 
an adiabatic temperature stratification and efficiency as given by Eq. [TSi 



Oscillation frequencies are computed using the Liege Oscillation code (LOSC. IScufiaire et al. 



20081 ) . The best-fit model is found via optimizations following the Levenberg-Marquardt 
algorith m, using spectroscopi c constraints and matching individual frequencies after apply- 
ing the iKjeldsen et al.l (120081 ) surface correction. The free parameters adjusted during the 
procedure are mass, age, initial helium abundance and metallicity, and the convectiv e and 
overshooting efficiencies. The same procedure was used by lMiglio fc MontalbanI (120051 ) using 
the frequency ratios ro2 instead of individual frequencies. 



5.1.5. Models E 



Models E were comp uted using the Asteroseismic Modeling Portal (AMP, e.g.. 



2012 



Metcalfe et al. 



Mathur et al.ll2012l ). It is based on a genetic algorithm for minimization ( Metcalfe &: Charbonneau 



2OO3I ) that matches the observed individual frequencies and spectroscopic constraints to those 



produced by a stellar evolution and pulsation codes. Its current implementation uses ASTEC 
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for the evoluti onary calculations, with simila r input physics as models B but low-temperatur e 
opacities from [Alexander fc FergusonI (119941 ) and nuclear reactions from lBahcall et al.l ( 119951 ) . 
No overshooting is included. 



Theoretical frequencies are computed with ADIPLS and corrected using the lKjeldsen et al, 

( I2OO8I ) formulation. The best-fit model is chosen by minimizing the sum of two separate x^ 



values: one determined for the individual frequencies and another one for the spectroscopic 
constraints. 



5.1.6. Models F 



Models F were computed using the Yale Rotation and Evolution Code ( JDemarque et al. 



20071). in its non- rot ating configuration. The code uses low-temperature opa cities from 

Ferguson et al.l ( l2005l ) and OPAL opacities for h igh te mperatures, OPAL EOS, the lGrevesse fc Sauval 
(119981 ) solar mixture, and lAdelberger et al.l (19981) n uclear reactions using the updated 

Forraicola et al.l (|2004l ). Models include helium and heavy- 
1994 ). and overshooting as given by Eq. [16] using a fixed 



14 



N(p, 7)^^0 cross section froni 



element diffusion (iThoul et al. 
value of Gov = 0.2 and an adiabatic temperature gradient when a convective core exists. 
This value of ov ershooting efficiency is m otivated by the fit to the open cluster NGC 2420, 
and described in iDemarque et al.l ( 119941 ) . 



The atmospheric stratification is given by an Eddington T — r relation, and frequencies 
calculated with the pulsation code used by lAntia fc Basul (119941 ). In order to find the best 
fit to the data, sequences of models were calculated with masses between 0.94-1.12 Mq for 
Perky and 1.08-1.24 Mq for Dushera. For each mass, we considered three different values 
of the mixing-length parameter, amit =1.55, 1.7, and 1.826 (the solar calibrated value). 
The best fit was found using the reduced y^ determined for the individual frequencies and 
frequency ratios as suming the solar surface term described in the supplementary material of 
Carter et all J2OI2I ). 



5.1.7. Models G 



Models G were constructed using ASTEC with the same input physics as models B, 
exploring different values of overshooting efficiency ctov = 0.0, 0.1, 0.2 and metallicity —0.5 < 
[Fe/H] < 0.5. No diffusion of elements was included in these calculations. A grid of models 
spanning a wide region of the Hertzsprung-Russell diagram was considered, with masses 
varying from 1.0 to 1.6 Mq. ADIPLS was used to compute the oscillation frequencies and 
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the best models were selected solely based on their asteroseismic characteristics, in partic- 
ular, their ro2 and tqio ratios, the frequency derivatives (slopes) of the ratios and the large 
frequency separations. In short, the range of possible models was first reduced based on the 
following criteria: that the slopes of ro2 and roio were within the l-a observed values and 
that the large frequency separation laid between the observed value ± 3 /iHz. The identifi- 
cation of the best model within this subset was based on the minimum of the x^ computed 
by comparing the ratios ro2 and roio with their observed values. 



5.1.8. Models H 



Mo dels H were extracted from a grid of models calculated using GARSTEC, described in 
detail in ISilva Aguirre et al.l ( 120121 ). Oscillation frequencies were computed using ADIPLS, 
and the best fits to the dat a were obtained looking for the minimum x^ of the individual fre- 
quencies after applying the lKjeldsen et al.l ( l2008l ) surface correction. No further optimization 
beyond the pre-computed grid was made. Input physics is the same as f or models A with th e 
following exceptions. The EOS is the FreeEOS ( A. W. Irwirg, see alsolCassisi et al.ll2003l ). 
hydrogen burning nuclear reaction rates are from lAdelberger et al.l (l201ll ). and amit = 1.811 
is fixed from the calibration of a standard solar model. Overshooting is implemented as in 
Eq. [15] with a fixed efficiency of ^ = 0.02 (also constrained geometrically in the case of small 
convective cores and thin shells) and microscopic diffusion of helium and metals has been 
included in all these models. 



5.2. Applying the surface correction 

In the previous sections we showed that search strategies and criteria for defining which 
model best reproduces the data are not unique and vary among teams. Some consider only 
individual frequencies while others also take into account frequency combinations. 

Figure [3] depicts the effect of applying the surface correction in one of the models of 
Perky (cf. Eq [T5]) . The upper panel shows the difference between corrected and uncorrected 
frequencies, presenting the same functional form as the solar case to which it is calibrated. 
The differences in the region of observed modes can be of up to ~ 20 /iHz, a factor of 100 
or more larger than the typical uncertainties in frequencies. 

A different picture arises when looking at the difference between the ratios tqio of cor- 



^http://freeeos. sourceforge.net/ 
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rected and uncorrected frequencies (lower panel in Fig. [3]). Discrepancies are at the 10" 
level, three orders of magnitude below the average uncertainties, s howing that this combi- 



natio n effectively cancels out the contribution from the outer layers (JRoxburgh fc Vorontsov 



20031 ). The value of the reduced x^i when considering only individual frequencies, is thus 



heavily affected by the use of the surface correction, while the same is not true for the ratios. 



Those teams applying the empirical surface correction of iKjeldsen et al.l ( l2008l ) consider 
somewhat different implementations of it. Before discussing how the parameters involved in 
the correction are actually obtained by each modeling team, a few words on the assumptions 
behind this formulation are necessary. It requires the value of r to be close to unity, meaning 
that the mean stellar density of the model actually reproduces that of the real star. It 
considers only radial modes when finding the parameters for the corrections, assuming that 
non-radial frequencies are equally affected by the surface layers. Finally, usage of a single 
power law is justified by considering the offsets in the mode frequencies to be caused solely 
by the properties of the outer stellar layers, and thus dependent on surface gravity, effective 
temperature, and surface composition. As a consequence of these assumptions models found 
using this empirical treatment are expected to reproduce the global stellar properties (mass 
and radius), but not necessarily the co re structure nor provide reliable age determinations 



see Section 5 in iKjeldsen et al.l l2008l ). Considering also non-radial modes in the fit is 



therefore necessary to ensure a proper representation of the deep stellar interior. 

Obtaining the parameters involved in the surface correction deserves some comments. In 
principle, for each combination of evolutionary and pulsation code the power-law exponent 
h should be determined from a solar calibration. This would also define the input physics 
to be used when modeling the targets, including the mixing-length p arameter. In practice , 



variations of h do not affect the obtained mean density of the models (IKjeldsen et al.ll2008l ). 
and one unique value can be used to model all stars. The h parameters applied by the 
modeling teams in this study have been obtained either via a solar calibration of their own 
codes or taken from others. Therefore the adopted values vary among the results. 

Another point where differences arise is in the calculation of the scaling parameter r. 
The best model for each star should have an r value very close to unity. In the original 
formulation, this is determined from theoretical and observed frequencies of radial modes 
only. Although it is the most commonly adopted manner, some modeling teams do not 
consider this r-scaling, or they estimate r using also the non-radial modes (models B in 
Section [5.1.21 above) . Finally, the reference frequency z/q is arbitrarily chosen. 

Since the surface correction is an empirical law based on the observed departure between 
the frequencies of the Sun and a solar model, there is no unique way of applying it. The 
net result is that x^ values for individual frequencies determined by each team are not 



-22- 



directly comparable due to these different implementations. Thus, we have uniformly applied 
the surface co r rectio n to all results to homogeneously compare them. As suggested by 



Kjeldsen et al.l (|2008[ ). we fixed b = 4.9, calculated r only with i = modes, and considered 
the observed frequency of maximum power Uuiax from each star as the reference frequency 
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Fig. 3. — Upper panel: frequency differences between corrected and uncorrected modes of 
oscillation. Typical uncertainties in the individual frequency determinations are of the order 
of ~0.2/xHz. Lower panel: differences between corrected and uncorrected frequency ratios 
roio- Typical uncertainties in the ratios are larger than the y-scale of the plot. See text for 
details. 
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5.3. Convective core location 

In Section H] we described liow tlie presence of a convective core produces an abrupt 
change in the adiabatic sound speed. The sharpness and shape of that transition will depend 
on the mixing processes taking place and the way overshooting is treated, if considered at 
alL In Fig. m some examples of these transition regions are shown for models with convective 
cores, where distinctive features can be observed in different cases. 

Model B in Fig. H] does not include overshooting and shows simply a sharp vertical 
transition at the edge of the core. Model G on the other hand applies a chemically homo- 
geneous extension of the convective region when applying overshoot, and thus presents a 
similar structure beyond the formal limits of the convective core. In contrast model F shows 
a smoother profile and a more extended transition region, produced by the interplay between 
overshooting and microscopic diffusion at the edge of the core. 

Since different assumptions and evolutionary codes produce different sound-speed pro- 
files, it is necessary to determine one criterion defining the position of the convective core for 
all models regardless of the considered input physics. As the p-modes propagate according 
to the sound speed, they are most sensitive to the regions where rapid variations in it take 
place. Therefore, we use the maximum in the adiabatic sound-speed derivative as a com- 
mon definition for the boundary of the homogeneously mixed central region. This basically 
defines the extent of the mixed zone, which is equivalent to the formal convective core size 
unless overshooting is included. 



6. Results 

In the following sections, we compare the best-fitting models found by each modeling 
team as described in Section 15.11 It is worth mentioning that, when comparing radius 
determinations, we decided to use the radius measured at the photosphere in order to avoid 
systematics introduced by the different prescription for atmospheric stratification used in 
each evolutionary code. 

In all cases, the values of xl for frequencies and Xr02 ^^^ Xroio fo^ frequency ratios 
(Eqs. [5] and |6]) given in the results below are calculated using frequencies corrected as ex- 
plained in Section 15.21 For comparison purposes we also give the results of the fit to the 
individual frequencies using the original frequencies as corrected by each modeling team, 
which we denote X^orig- Since the frequency correction does not affect the computed ratios, 
we will sort models according to the average between Xroio ^^^ Xr02- We remind the reader 
that the correlations in the ratios roio are taken into account when computing XroiO' ^^ 
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Fig. 4. — Adiabatic sound speed near the stellar core as a function of acoustic radius for a 
set of selected models: B (red solid line), F (violet dotted line) , and G (black dashed line). 
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mentioned in Section [5l 



6.1. Perky 



Table 3: Main characteristics of the models calculated for Perky from each method as de- 



)ed in Section 15.11 


















Model 


M{Mq) 


Age (Gyr) 


logs 


Teff (K) 


Rphot/RQ 


[Fe/H^ 


z^ 


Y^ 


log (L/Lq) 


Model A 


1.09 


4.62 


4.299 


6074 


1.2257 


-0.09 


0.0181 


0.273 


0.263 


Model B 


1.11 


4.36 


4.305 


6122 


1.2284 


+0.00 


0.0175 


0.267 


0.279 


Model C 


1.05 


5.53 


4.294 


5970 


1.2093 


-0.09 


0.0145 


0.260 


0.222 


Model D 


1.17 


4.97 


4.306 


6016 


1.2570 


-0.01 


0.0210 


0.240 


0.268 


Model E 


1.17 


4.94 


4.307 


5951 


1.2580 


-t-O.Ol 


0.0197 


0.226 


0.250 


Model F 


1.12 


4.88 


4.301 


6100 


1.2370 


-0.09 


0.0179 


0.264 


0.279 


Model G 


1.06 


4.47 


4.299 


6071 


1.2088 


-0.10 


0.0139 


0.271 


0.250 


Model H 


1.13 


3.80 


4.305 


6134 


1.2382 


-0.07 


0.0188 


0.269 


0.289 



Table 4: Structure and seismic characteristics of the models for Perky. Values of x^orig were 



obtained using the frequencies as corrected by each modeling team, while xi was calculated 
using the standard version of the surface correction. Last column gives the symbol used to 
plot each model in all figures. 



Model 



Xc 



ice (s) iphot (s) tTot (s) tcc/tphot (%) xSc 



v2 

ArOlO 



v2 
Ar02 



Symbol 



Model A 


0.124 


Model B 


0.249 


Model C 


0.078 


Model D 


0.085 


Model E 


0.120 


Model F 


0.085 


Model G 


0.122 


Model H 


0.195 



95.0 



4549.5 


4577.4 


4541.9 


4621.7 


4557.5 


4593.7 


4565.3 


4691.6 


4573.0 


4633.3 


4553.4 


4725.7 


4546.3 


4626.1 


4547.6 


4613.6 



2.09 



- 


21.28 


0.79 


0.35 


Filled circle 


9.15 


8.96 


1.07 


0.43 


Open square 


11.03 


15.40 


1.52 


1.12 


Filled cross 


2.37 


3.25 


1.67 


1.21 


Filled square 


5.04 


6.05 


2.23 


1.39 


Downward triangle 


3.50 


10.27 


1.04 


3.62 


Open diamond 


- 


21.78 


2.22 


3.56 


Upward triangle 


19.19 


20.13 


1.97 


5.46 


Filled diamond 



The main characteristics of the best-fit model found by each team for the target Perky 
are given in Table [31 The first feature to notice is that the agreement in surface gravity 
between all models is better than 0.02 dex. This level of precision is already better than 
what was found using the direct and grid-based methods (see Table H]), and the obtained 
mass range from detailed modeling is now restricted to values between 1.05-1.17 M©. This 
decrease is partly responsible for the reduction in the age difference compared to the grid 
results, now ranging from ~3.8 to ~5.5 Gyr. The age spread is constrained within a level 
of ~18%, which although not extremely good is much better than that obtained from the 
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global asteroseismic fit (see Section [3]). Effective temperatures are also contained within the 
values given in Table [H The overall agreement in radius is of the order of 0.05 Rq, only 
slightly better than that found with the global analysis. 

Models D and E have initial helium abundances below the currently accepted SBBN 
value. This problem occurs regularly when models are selected solely on the base of best 
reproducing the observed individual freq uencies after applying the surface correction (e.g., 



Metcalfe et al.ll2010l . iMathur et al.l 120121 ). Interestingly enough, models D and E have basi- 
cally the same characteristics. Both cases present the highest mass values among the results, 
hinting that the correlation between mass and helium abundanc e through the mean m olec- 



ular weight is reflected in the asteroseismic fits (see also Fig. 5 in lMetcalfe et al.ll2009[ ). The 
initial metallicities of these models are also the highest ones, probably to compensate for the 
low helium abundance when achieving a similar luminosity and effective temperature. This 
might point towards deficits in the use of opacity or EOS data. We discuss this further in 
Section [71 

As we aim to make a better characterization of the target properties we analyze the 
spectrum of oscillations in detail. In Table H] the internal structure and asteroseismic prop- 
erties are given for the best-fit models of Perky. The values of Xroio have been calculated 
using the frequency range marked in Fig. [H Interestingly, the models with the lowest xl 
value for the individual frequencies are not the ones that best reproduce the observed ra- 
tios, and are in fact those with the lower than SBBN initial helium abundances (D and E). 
Model F presents a large difference between the original surface correction x^orig ^^^ the one 
we applied (y?,). T^i^ rn odel was the only one selected using a prescription different from 



the iKjeldsen et al.l ( l2008l ) one (see Section [5^ . reflecting the large impact of the assumed 
recipe for surface correction and warning us about selecting models solely based on fits to 
individual frequencies. 

A more detailed picture can be made when looking at Figures [5] and O where we show 
the comparison to the data of models A, E, F, G and H. From Fig. [5] it can be seen that the 
ratios tqio of model E (G) are systematically below (above) the observed ones, allowing us 
to discriminate against them. Ratios ro2 in Fig. IH] have larger error bars and less data, but 
those of model G are also in disagreement with the observations. Following the same line 
of reasoning models F and H can also be discarded due to their large discrepancy with the 
observed values of ro2. 

More interesting than discarding models is trying to understand the structural differ- 
ences that produce this behavior. As discussed in Section HJ the ratios ro2 are a proxy for 
age and thus central hydrogen content of the star. Inspection of Table H] does not reveal a 
correlation between Xc and the goodness of fit to ro2. In fact, models with the same central 
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hydrogen content (D and F) present very different values of Xr025 wliile model B gives an 
excellent fit to these ratios with more than twice the amount of hydrogen remaining in the 
core than any other model. Because of the systematic differences in the input physics and 
evolutionary codes, we no longer find a correlation between ro2 and Xc- 

The case of model B presents an interesting challenge to our current understanding of 
asteroseismic diagnosis. With the exception of this one, all other results point towards a 
star w ith no convective core . Nonetheless, this model was selected using a modified version 



of the iKjeldsen et al.l (120081 ) surface correction that includes also the non-radial modes (see 
Section |5^ . naturally resulting in a xl value that is slightly different than what would be 
obtained using the radial modes only. Interestingly enough, the algorithm applied to select 
model B returns a different best-fit case when the "normal" surface correction is applied 
(using only i = modes). This new solution, which we call model B', gives values of 
xloTig = 9-27 and xl = 8.18. When compared to the results of model B given in Table H] 
of xloTig = 9-15 and xl = 8.96 we see that model B' was not selected because it used the 
modified surface correction, while it would have been chosen if the standard one was applied 
instead. 

Model B', has exactly the same initial composition as model B in Table |3l but with 
a mass of 1.10 Mq and no convective core present. Its age is 4.72 Gyr, Tcs =6105 K, 
log g = 4.300 and luminosity log (L/Lq) = 0.275, in good agreement with the results of 
model A and similar to model B. However its central hydrogen abundance of Xc = 0.03 
is much lower than found in the other solutions, and gives fits to the frequency ratios of 
Xroio = 3-55 and xlm = 0-43. 

Figure [7] shows the comparison to the data of models A, B, and B'. Despite the very 
different central hydrogen contents (0.115, 0.249 and 0.03 for models A, B, and B', respec- 
tively), the ratios ro2 show no significant sensitivity to this parameter. Inspection of the 
ratios roio allows us to discriminate against model B', but it is clear that for this particular 
case our current diagnostic tools are not capable of discriminating against the presence of 
a convective core (i.e. between models A and B). We remind the reader that using the 
standard formulation of the surface correction from Eq. [T^ would have resulted in model B' 
as the best fit found by this team, a model clearly not compatible with the observed ratios 
'^oio- 

To explore this issue in more detail, we reproduced the results of model B with the 
evolutionary code and input physics of model A. For the same initial metallicity and mass 
it was possible to sustain a convective core beyond the pre-main sequence phase only by 
removing the geometrical restriction for small convective cores (see Section I5.1.ip . The 
resulting model has an age of 4.71 Gyr at Tefj = 6148 K, log g = 4.303 and luminosity 
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log [L/Lq) = 0.288. Convective core size is slightly smaller than the one found in model B, 
placed at an acoustic radius of t = 81.2 s that corresponds to 1.77% of the total photospheric 
acoustic radius. The fit to the frequency data is also comparably good as in model B, with 
X?oio = 1-45 and Xr02 = 0-32. 

The overall picture we are left with is the following. Based on the analysis described in 
the previous paragraphs, we can discard models E, F, G, and H and obtain from the rest of 
the results a mean and standard deviation pointing towards a star oi R = 1.23 ± 0.02 Rq, 
M = 1.11 ± 0.05 Mq, and 4.87 ± 0.50 Gyr. This corresponds to uncertainties of 1.6%, 4.5% 
and 10.4% in radius, mass, and age, respectively; the level of precision obtained in the stellar 
parameters by this method of asteroseismic analysis is extremely good, particularly in age. 
However the experience of model B shows that slightly different models with equally good 
or even better fits to the data can still be lurking among our results, and have not been 
picked out due to our assumptions on the fitting algorithms. Although the mass, radius and 
age are not particularly affected by this issue, there could be other sources of systematic 
errors arising from unexplored input physics. In any case, we expect the true uncertainties 
to be only sightly larger given the heterogeneous set of codes and input physics used in this 
analysis. 

We are left with the unsolved issue of whether a convective core exists or not and, in 
relation to that, what the amount of hydrogen still remaining for central burning is. This 
has a direct impact in the inferred age at the turn-off: while model A will remain in the 
main sequence for ~1 Gyr, model B will spend another ~2 Gyr in this evolutionary phase. 



6.2. Dushera 

For our second target of interest, Dushera, we performed a similar analysis to the one 
made for Perky. In Table O we give the main characteristics of the best-fit model found by 
each team, ranked again by their average fit to the frequency ratios. In this case no model 
gives a result with an initial helium abundance below the SBBN value. Agreement in log g 
values is again better than 0.02 dex, while the radii of the results are contained within a 
range of 0.06 Rq. 

The mass range of the models is similar to the one obtained via the grid-based method, 
but with an age spread of the order of ~ 15%. All analyses for this target returned models 
where a convective core exists, but there are differences in the size of the centrally mixed 
region that we aim to characterize using the frequency combinations. Table IH] presents the 
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Fig. 5. — Frequency ratios roio('T-) for a selected set of models for Perky, and data only in the 
range marked in Fig. [TJ Observational frequencies of Perky are depicted in open circles, while 
models A, E, F, G and H are shown as filled circles, downward triangles, open diamonds, 
upwards triangles, and filled diamonds, respectively. 



internal structure characteristics of each best-fit model, as well as the values of the fits to 
individual frequencies and ratios. 

Once again, inspection of Table [6] reveals that models best fitting the individual fre- 
quencies after surface correction (D and B) do not necessarily give the best match to the 
frequency combinations. Just as in the case of Perky, the data warn us to be cautious when 
relying in results solely based on fits to individual frequencies. Also, comparison of the values 
given in Table [6] shows no significant correlation between the central hydrogen content Xc 
and ro2 among these models built with different input physics. 



-30- 



0.080 



0.075 



0.070 



o 



Perky 




0.065 

•Model A 

^ Model E 

Model F 

0.060- A Model G 

♦Model H 

I I I I I I I 



1800 1900 2000 2100 

u (/i.Hz) 



2200 



2300 2400 



Fig. 6. — Frequency ratios ro2 for the same models of Perky as in Fig. [5l Data for only those 
modes considered in the fit are plotted (see Sect. S]). 



We can further constrain the properties of the target by analyzing the fits to the ratios, 
plotted in Fig. [8] for the three models presenting the largest average of Xroio ^^^ Xr02- 
Models H is systematically below the observed ratios roio, while model B shows a slope that 
is not consistent with the data. Moreover, when looking at the ratios ro2, all models are 
systematically lower than the observed values. 

In Fig. |9] we plot the ratios roio of the first 5 models in Table [51 The slopes of mod- 
els C and E are clearly different from the observed one, with model C (E) showing lower 
(higher) values than the data in the low (high) frequency range of the considered ratios. 
Model D presents a similar behavior as model C 
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Fig. 7. — Upper panel: frequency ratios ro2 of Perky. Lower panel: frequency ratios roio('^) 
of Perky. Observational frequencies are depicted by open circles, while models A, B, and B' 
are shown as filled circles, open squares, and open crosses, respectively. 



Models A and F are the only two that consistently fit all the available data of Dushera, 
reproducing also the slope of the observed ratios tqio and ro2. Considering only these results, 
our analysis points to a star with R = 1.39±0.01 Rq, M = 1.15± O.O4M0, and 3.80±0.37 Gyr 
corresponding to uncertainties of 1.0%, 3.1% and 9.9% in radius, mass and age, respectively. 
As in the case of Perky we caution the reader that these values are to be taken as a reference, 
since systematics can arise from unexplored input physics and from models with slightly 
different properties that can also fit the data correctly but are overlooked by the fitting 
algorithms. We have shown that models fitting all asteroseismic combinations are restricted 
to a very small stellar parameter space and thus we expect the total uncertainties to be only 
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Table 5: Main model characteristics of the best fits to the data of Dushera. Models are sorted 
according to the value of their Xroio ^^ given in Table |6l Last column gives the value of the 
overshooting efficiency parameter used. 



Model 


M (Mq) 


Age (Gyr) 


logs 


Toff (K) 


J2phot/R0 


[Fc/H^ 


z^ 


Y^ 


OV 


Model F 


1.12 


3.53 


4.208 


6318 


1.3803 


-0.15 


0.0162 


0.310 


Qov = 0.2 


Model A 


1.17 


4.06 


4.215 


6305 


1.3990 


-0.12 


0.0135 


0.261 


i = 0.03 


Model D 


1.26 


3.09 


4.223 


6305 


1.4378 


-0.03 


0.0202 


0.266 


Qov = 0.015 


Model E 


1.19 


3.44 


4.217 


6262 


1.4083 


-0.11 


0.0147 


0.263 


«ov = 


Model C 


1.20 


3.34 


4.218 


6094 


1.4119 


-f0.09 


0.0210 


0.280 


CKOV = 


Model G 


1.14 


4.28 


4.208 


6040 


1.3926 


-1-0.00 


0.0173 


0.278 


Qov = 0.1 


Model B 


1.20 


3.67 


4.216 


6205 


1.4157 


-1-0.00 


0.0175 


0.267 


Cfov = 


Model H 


1.19 


4.00 


4.214 


6150 


1.4119 


-0.08 


0.0188 


0.269 


i = 0.02 



Table 6: Internal structure and seismic characteristics for the models of Dushera. Values of 
X^orig were obtained using the frequencies as corrected by each modeling team, while x^ was 
calculated using the standard version of the surface correction. 



Model 


X, 


tec (s) 


Afcc (%) 


i?cc (%) 


tphot (s) 


tTot (s) 


tcc/tphot (%) 


X^ ■ 

'^v orig 


xl 


XrOlO 


Xrt2 


Model F 


0.236 


129.5 


8.32 


7.10 


5367.1 


5577.5 


2.41 


7.09 


11.81 


0.92 


0.93 


Model A 


0.262 


131.4 


7.94 


7.08 


5368.1 


5401.2 


2.45 


- 


20.3 


1.44 


0.47 


Model D 


0.256 


125.1 


5.89 


6.57 


5387.4 


5546.1 


2.32 


2.04 


1.90 


1.60 


0.53 


Model E 


0.213 


110.1 


4.77 


5.79 


5276.0 


5347.2 


2.09 


4.78 


4.67 


1.74 


0.72 


Model C 


0.236 


116.0 


5.01 


6.09 


5390.8 


5433.9 


2.15 


9.93 


10.31 


2.06 


0.64 


Model G 


0.214 


126.2 


6.94 


6.63 


5396.9 


5498.1 


2.34 


- 


9.16 


0.81 


2.91 


Model B 


0.173 


117.3 


5.28 


5.99 


5319.5 


5417.1 


2.21 


4.72 


3.87 


1.58 


2.30 


Model H 


0.243 


122.0 


6.66 


6.67 


5373.7 


5430.9 


2.27 


20.98 


16.93 


1.92 


4.10 



slightly larger than those given above. 

The analysis of Dushera not only reveals the presence of a convective core but it also puts 
constraints on its size. Our diagnostics are sensitive to the total extent of the mixed region, 
which in this case was reproduced by inclusion of overshooting. It is clear from Table [6] and 
Fig. [9] that only models with a mixed core comprising more than 2.35% of the total acoustic 
radius are able to reproduce the ratios. The overall results leave us with a homogeneously 
mixed region of ~2.4% t/tphot, and in any case extending beyond the formal convective 
boundary due to the inclusion of overshooting. Thus, we measure here for the first time 
mixing beyond the formal convective boundaries in the core of a Kepler field main-sequence 
star. 

Tests were made on the impact of the core size by running model A with higher and 
lower overshooting efficiency values to produce different sizes of mixed zones. These models 
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were not able properly to reproduce the ratios roio and show that, unlike the case of Perky, 
these ratios are mostly sensitive to the size of a convective core when this is already well 
developed during the main-sequence phase. We also performed tests using the model F set 
of evolutionary calculations on the impact of reducing the considered number of ratios tqio 
by two, and found no significant difference in the obtained solutions. 
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Fig. 9. — Ratios roio of Dushera. Selected models plotted are A (filled circles), C (filled 
crosses), D (filled squares), E (filled downward triangles), and F (open diamonds). 

7. Discussion 



Asteroseismic analysis aims to provide accurate stellar parameters of pulsating stars by 
theoretically reproducing their frequency spectrum. This has usually been done by match- 
ing the observed oscillation frequencies to theoretically calculated ones after applying an 
empirical correction to account for our poor modeling of stellar outer layers. Our results 
confirm that best-fit models found by this method do not always reproduce frequency com- 
binations, and that different criteria to define what the best-fit model is can lead to solutions 
with similar stellar pa rameters but very different interior structures. As originally stated by 
Kjeldsen et al.l ( l2008l ). the surface correction only ensures an accurate reproduction of the 
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stellar mean density by assuming that the frequency offsets are caused solely by the proper- 
ties of the outer layers. Usage of non-radial modes and frequency combinations is necessary 
to match properly the deep interior and provide a reliable estimate of stellar age. 

Of particular importance is the implementation of the surface correction itself, which has 
a direct impact in the goodness of fit via the reduced xl value. The values of the corrections 
are much larger than the errors in the frequencies, making the definition of the best-fit 
model sensitive to the way these corrections are applied. Thus, the xl criterion cannot 
be objectively used among best-fits found using different implementations of the surface 
correction. We have proposed instead to characterize the goodness of fit to asteroseismic 
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data by using the frequency ratios roio and ro2. These combinations are thought to be mostly 
determined by the central conditions of the star and almost unaffected by the outer layers, 
providing a more robust approach than matching individual frequencies. The masses and 
radii obtained by both techniques agree within 2-o", but the resulting models can have very 
different structural properties. Thus, although the overall results are not hugely affected by 
the fitting techniques and the issue of the surface correction, using frequency ratios instead 
of just reproducing the individual frequencies of oscillation is a more rehable way to yield 



accur ate stellar parameters, in particular age (as was also suggested by iMiglio fc Montalban 



20051 ). 



The frequency ratios have been used as proxies for different internal properties of stars. 
For instance, a tight correlation between ro2 and central hydrogen content (age) holds for a 
given set of input physics, as shown by the modified C-D diagram in Fig. [2J However, no 
correlation is found between these parameters among the set of models that fit the seismic 
data. The reason for this is that the equivalence between central hydrogen content and ro2 
i s only valid for a fixed set of input phy sics, which is naturally not the case in this study 



(IMonteiro et al.lbOOl IWhite et all 120111 1. 



In a similar way the ratios roio, which were thought to be sensitive to the existence 
and size of the convective core as well as a proxy of the central hydrogen abundance, were 
found to be affected by other parameters not yet fully identified. The latter sensitivity 
is particularly important in the case of Perky, where analysis of tqio did not allow us to 
discriminate between models with or without a convective core. The stellar parameters 
obtained at this stage are not significantly affected by the presence of a convective core; thus 
we can still derive a reliable mass, radius and age for this star. Nevertheless the remaining 
central hydrogen content of the target is not well constrained, and therefore neither is the 
age it will have at the end of the main sequence. 

More importantly, the experience of Perky is a warning flag to our fitting algorithms 
and the use of the surface correction. It is easy to imagine that if all teams used a different 
implementation of the surface correction, different solutions from the ones presented here 
would have been found. In fact, the restriction to the amount of overshooting in the case 
of small convective cores might be overestimated in cases like Perky. It has been shown for 
the solar case that restricting t he amount of pre-main sequence overshooting is necessary 



to reproduce helioseismic data ( ISchlattl fc Weisslll999l ). but for stars slightly more massive 



than the Sun this needs to be further explored. 

Related to fitting techniques, it remains fully to understand why in some cases low initial 
helium abundances are obtained when only individual frequencies are matched. Although 
the results for Perky suggest that mass is the major contributor to this feature, the overall 
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picture is certainly more complex. While there is no reason to believe a universal value of 
AY/ AZ applies for all stars, and small variations in helium are therefore expected, sub- 
SBBN values are clearly a challenge to stellar models. Such a be havior already appeare d 



from the analysis of the main-sequence broadening in field K dwarfs (ICasagrande et al. 



20071 ) 



and sim ple analytical conside rations suggested that model opacities could have been a likely 
culprit ( jPortinari et al.ll2010l ). 



Mathur et al.l ( 12012| ) performed a uniform asteroseismic analysis of 22 Kepler targets us- 



ing AMP, the same code and fitting algorithm of models E described in Section [5.1.51 Their 
final results show that 10 of their 22 best-fit models have a mean value of initial helium 



below that of SBBN (see their Table 5). The models in lMathur et al.l ( 120121 ) are not directly 
comparable to ours since they are based on analysis of only 1-month long observations, but 
the overall outcome clearly points out to the need of further studies in this topic. These are 
currently underway (J0rgensen et al. 2013, in preparation). Understanding these discrep- 
ancies is important not only for the sake of stellar models, but for the possibility of using 
these stars as tracers of the Galactic AY/ AZ and possibly one day for multiple populations 
in stellar clusters 

Stars like Perky lie at the limiting mass where the onset of core convection is thought 
to occur. A lthough such targets ha ve been used in open clusters as a test for the solar 
abundances (IVandenberg et al.l 120071 ). it has been shown that the existence or absence of 
a convective core in stellar models heavily depends on our assumptions about the input 
physics. The upcoming releases of Kepler data will give us more frequencies with higher 
accuracy and we can certainly hope to solve this issue. 

In the case of Dushera, we have made the first direct detection of a convective core in a 
iTep/er main-sequence target. Use of fits to the ratios roio and ro2 allow us to discard outliers 
and determine the characteristics of the star with an excellent accuracy. The size of the 
central mixed region is estimated to be ~2.4% of the photospheric acoustic radius, showing 
that mixing beyond the formal Schwarzschild convective boundary exists. In this study 
the extra mixing has been modeled using different overshooting prescriptions; in reality it 
could also be produced, for instance, by rotational mixing or a combination of more physical 
processes. 

For both targets, masses and radii obtained by our asteroseismic analysis have a preci- 
sion of ~4% and ~1.5%, respectively. Regarding stellar ages, we determine them with an 
encouraging ~10% level of precision. This is probably optimistic due to our fitting algo- 
rithms potentially overlooking models that might better reproduce the data, and also since 
our models do not explore all possible sets of input physics. Thus, additional sources of 
uncertainties can exist that have not been explored yet. In fact, systematics arising from. 
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for instance, the unknown input p hysics can easily account for an extra ~7% in age at the 
end of the main-sequence phase (IValle et al.l l2013l ). Although our analysis confirms that 



we are certainly doing better than the ~20-50% in age obtained with the usual isochrori e 



fitting from a single stellar evolution database (e.g., iTakeda et al.l 120071 . ISoderblomI |2010| ) . 
asteroseismology has the potential to deliver even better ages once we fully understand our 
current diagnostic tools and pursue further theoretical studies to develop even better ones. 
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Table 7: Frequencies of Perky 



i 


i/ (pViz) 





1705.262 ±0.150 





1807.695 ±0.195 





1909.986 ±0.151 





2012.862 ±0.132 





2117.342 ±0.113 





2221.532 ±0.103 





2325.586 ±0.136 





2429.842 ±0.151 





2533.811 ±0.206 





2639.023 ±0.351 





2744.557 ± 0.778 
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1957.109 ±0.125 




2061.479 ±0.107 




2165.840 ±0.110 




2270.400 ±0.116 




2374.555 ±0.139 




2479.027 ±0.178 




2583.464 ± 0.213 




2688.717 ±0.362 




2794.513 ±0.515 


2 


1697.972 ± 0.254 


2 


1800.063 ± 0.392 


2 


1902.316 ±0.506 


2 


2005.719 ±0.184 


2 


2110.165 ±0.114 


2 


2214.715 ±0.202 


2 


2318.862 ±0.180 


2 


2422.900 ± 0.288 


2 


2528.083 ± 0.373 


2 


2631.927 ±0.700 


2 


2740.448 ± 0.890 



